1 
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The local void model has recently attracted considerable attention because it can explain 
the apparent accelerated expansion of the present universe without introducing dark energy. 
However, in order to justify this model as an alternative to the standard /1CDM cosmology, 
the model should be tested by various observations, such as the CMB temperature anisotropy, 
besides the distance-redshift relation of SNIa. For this purpose, we derive analytic formulae 
for the dipole and quadrupole moments of the CMB temperature anisotropy that hold for 
any spherically symmetric universe model and can be used to compare consequences of such 
a model with observations of the CMB temperature anisotropy rigorously. We check that our 
formulae are consistent with the numerical studies previously made for the CMB temperature 
anisotropy in the void model. We also update the constraints concerning the location of the 
observers in the void model by applying our analytic dipole formula with the latest WMAP 
data. 

§ 1 . Int r o duct ion 

In modern cosmology, it is commonly assumed that our universe be isotropic 
and homogeneous on large scales and accordingly be described by the Friedmann- 
Lemaitre- Robertson- Walker (FLRW) metric as a first approximation. Then, together 
with consequences of various cosmological observations such as the power spectrum of 
Cosmic Microwave Background (CMB) temperature anisotropy, the distance-redshift 
relation of type la supernovae (SNIa) indicates that the expansion of the present uni- 
verse is accelerated. We are then led to introduce "dark energy," which has negative 
pressure and behaves just like a positive cosmological constant, thus having the stan- 
dard ylCDM model. However, there does not appear to be any satisfactory theory 
that can naturally explain the origin of dark energy and its magnitude required by 
observations. It is therefore tempting to seek alternative explanations of the ap- 
parent cosmic acceleration, more specifically, the SNIa luminosity distance-redshift 
relation, without invoking dark energy. (For a number of approaches to the dark 
energy problem, see e.g., 1) and references therein). 

One of such attempts is a "local void model" proposed by Tomita, 2 )' 3 ) also in- 
dependently by Celerier, 4 ) and by Goodwin et al. 5 ** In this model, our universe is 
no longer assumed to be homogeneous, having instead an under-dense local void in 
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the surrounding overdense universe. *•* The isotropic nature of cosmological observa- 
tions on large scales is realized by assuming the spherical symmetry and demanding 
that we live close to the center of the void. Furthermore, the model is supposed 
to contain only ordinary dust like cosmic matter, describing, say, the CDM com- 
ponent. Since such a spacetime can be described by the Lemaitre-Tolman-Bondi 
(LTB) metric, 23 )~ 25 ) we also call this model the "LTB cosmological model." **) Since 
the rate of expansion in the void region is larger than that in the outer overdense 
region, this model can account for the observed dimming of SNIa luminosity. In 
fact, recent numerical analyses 28 ) -32 )' 34 ) -40 )' 42 )' 43 ) have shown that the LTB model 
can accurately reproduce the SNIa distance-redshift relation. For this reason, de- 
spite the relinquishment of the widely accepted Copernican/cosmological principle, 
the local void/LTB model has recently attracted considerable attention. 

However, in order to justify the LTB model as a viable alternative to the stan- 
dard ylCDM model, one has to test this model by various observations other than 
the SNIa distance-redshift relation. A number of papers for this purpose have ap- 
peared recently, studying constraints from observations, such as the CMB tem- 
perature anisotropy, 31 ) -33 )' 43 ) the baryon acoustic oscillation, 41 )~ 43 ) the kinematic 
Sunyaev-Zeldovich effect, 44 )' 45 ) the galaxy count-redshift relation, 46 ) -49 ) etc. Many 
LTB models with a small, a few hundreds Mpc size void are already ruled out, but at 
present, the models with a huge, Gpc size void still remain to be tested. More details 
of the current status of the LTB models are discussed in 6). Most of these analyses 
have been performed for various types of LTB models by using numerical methods, 
and it does not seem to be straightforward to compare analyses for each different 
model so as to have a coherent understanding of the results. In order to have general 
consequences of the LTB cosmology and systematically examine its viability, it is 
desirable to develop some general, analytic methods that can apply, independently 
of the details of each specific model. 

Apart from the quest of alternatives to dark energy, the LTB metric may also 
serve as a toy model for getting some insights into the dynamics and possible observa- 
tional effects of non-linear perturbations in the standard ylCDM cosmology. The LTB 
metric can incorporate a cosmological constant in a straightforward manner, hence 
being able to describe, as an exact solution, highly non-linear inhomogeneities — 
almost arbitrary in magnitude, as long as being spherically symmetric — in the FLRW 
universe with dark energy. In view of this, it is also worth attempting to derive some 
analytic formulae that can be used to make theoretical predictions of the yl-LTB 
spacetimes and rigorous comparison with cosmological observations. 

The purpose of this paper is to derive analytic formulae for the dipole and 
quadrupole moments of the CMB temperature anisotropy in general spherically 



*' See, e.g., 6) for earlier work along a similar line. For other proposals that also attempt to 
exploit the effects of inhomogeneities to solve the dark energy problem, interests and criticisms 
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**' In fact, the LTB metric has been considered in the cosmological context even much before 
the present form of the dark energy problem was raised by the results of SNIa observations. For 
example, an early proposal of the void model constructed by the LTB metric was made by Moffat 
andTatarski. 26 )- 27 ' 
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symmetric inhomogeneous spacetimes, including the yl-LTB spacetime as a partic- 
ular case. In the standard cosmology, basic properties of the CMB temperature 
anisotropy are derived by inspecting perturbations of the Einstein and Boltzmann 
equations in the FLRW background universe. Ideally, it is desirable to do the same 
thing in the LTB background. However, perturbations in a LTB spacetime, let alone 
those in a general spherically symmetric spacetime, have not been very well studied 
mainly because the perturbation equations are much more involved to solve in the 
spherical but inhomogeneous background, though the linear perturbation formulae 
themselves have long been available 50 )' 51 )*). In this paper, we are not going to deal 
with perturbations of the LTB metric. Instead, we will exploit the key requirement 
of the LTB cosmology that we, observers, are restricted to be around very near the 
center of the spherical symmetry: Namely, we first note that the small distance be- 
tween the symmetry center and an off-center observer gives rise to a corresponding 
deviation in the photon distribution function. Then, by taking 'Taylor-expansions' 
of the photon distribution function at the center with respect to the deviation, we 
can read off the CMB temperature anisotropy caused by the deviation in the photon 
distribution function. By doing so, we can, in principle, construct the Z-th order 
multiple moment of the CMB temperature anisotropy from the (up to) l-th order 
expansion coefficients, with the help of the background null geodesic equations and 
the Boltzmann equation. We will do so for the first and second-order expansions to 
find the CMB dipole and quadrupole moments. We also provide the concrete expres- 
sion of the corresponding formulae for the LTB cosmological model. Our formulae 
are then checked to be consistent with the numerical analyses of the CMB temper- 
ature anisotropy in the LTB model, previously made by Alnes and Amarzguioui. 53 ) 
We then apply our formulae to place the constraint on the distance between an ob- 
server and the symmetry center of the void, by using the latest Wilkinson Microwave 
Anisotropy Probe (WMAP) data, thereby updating the results of the previous anal- 
yses. 

In the next section, we derive analytic formulae for the CMB temperature 
anisotropy in the most general spherically symmetric spacetime. In §3, we obtain 
analytic formulae for the CMB temperature anisotropy in the LTB model, and some 
constraints concerning the position of the observer. §4 is devoted to summary. Ap- 
pendix discusses the regularity of the photon distribution function and the behavior 
of some geometric quantities at the symmetry center. 

§2. The CMB temperature anisotropy in spherically symmetric 

spacetimes 

In this section, we will derive analytic formulae for the dipole and quadrupole mo- 
ments of the CMB temperature anisotropy. We first briefly discuss the null geodesic 
equations, the photon distribution function and its relation to the CMB photon tem- 
perature in the most general spherically symmetric spacetime. Next, inspecting the 
first-order derivatives of the photon distribution function at the center of the spher- 
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ical symmetry and using the null geodesic equations, we derive the dipole formula 
for the CMB. Then, taking further derivatives, we derive the quadrupole formula of 
the CMB temperature anisotropy in general spherically symmetric spacetimes. 

2.1. The photon distribution function in a general spherically symmetric spacetime 

The most general spherically symmetric metric can be written in the following 
form 

ds 2 = -N 2 (t,r)dt 2 + S 2 (t,r)dr 2 + R 2 (t,r)(d6 2 + sin 2 Odcp 2 ). (2-1) 
In this background, the relevant geodesic equations are given by 

d A - A {v t ? _ _ ss 2 _ _r_ , 

d\~ N {P) N PP N* [p) N*R P± > {22) 

^ - - — (J) 2 - - ~(v r ) 2 + —v 2 (2-3) 

dX ~ S 2 {P} S PP S {P} + S 2 R P ^ (23) 

where p^ = dx^/dX, p^p^ 1 = 0, with A being an affine parameter, and p± is given 
by p\ = R 2 {(p 9 ) 2 + {p^) 2 sin 2 6}. Here and in the rest of the paper, prime and dot 
denote the derivatives with respect to r and t, respectively. 

We are concerned with the distribution function, F(x,p), for the CMB photons 
that leave the "last scattering surface" appropriately defined, say, t = U hypersur- 
face, in the universe modelled by the above metric and that eventually reach an 
"observer" very near the symmetry center. We assume that the photon distribution 
function F(x, p) is spherically symmetric, respecting the symmetry of the background 
geometry, so that 

F{x,p) = F Q (t,r,u,n), (2-4) 

where uj = p l and n = Sp r /(Nuj). Note that from the above geodesic equations, we 
have 

N , ^ 

r = M^, (2-5) 

(N N' 9 S , n x i?\ , N 

U = ~ U \N +2fi ~S +tl S + { 1 ~ fi) R ' (2 ' 6) 



A = (l-/x 2 ) -r --T7 o - o ■ (2-7) 




In particular, it follows from the above equations 

du>r = dull = 0, d r fi = 0, (2-8) 

where the last one holds for the radial null geodesies, for which fi = ±1. 

Now, suppose the universe is locally in thermal equilibrium, that is, F is given 
as the Planck distribution function, 0, at the last scattering surface, as required in 
most of the known LTB models. Then, given a photon geodesic 7, the ratio of the 



5 



temperature T of the photon and its energy oj is preserved along the trajectory 7. 
Therefore, ui comes in F in the form 

F = <P(u/T). (2-9) 

Then, the CMB temperature anisotropy ST/T is generally given by 

{SF) W + {§F )W +... = j-^ + \ (f ) 2 M,) 2 + ■ ■ ■ J *. (2-10) 

Now, suppose that an observer lives at a distance of 5x l from the symmetry center. 
Then, the left-side are written as 

(8F)W = Sx'idiF^, (6F)W = ^Sx'SxHdidjFh, (2-11) 

where here and in the following, the subscript '0' implies the value evaluated at the 
center (r = 0) at the present time (t = to). The CMB temperature anisotropy dipole 
(8T/T)W and quadrupole (5T/T)W are therefore given by 



5T\ {1) _ S^jdjFjo 
~f) " cod^Fo ■ 



(2-12) 



6T\V> _ l Sx'SxijdjdjFh 1 \J5T\W\ {ud u )*F 
T J 2 ud u F ' 2\\T J j ud^Fo ' ( ' ' 

Since (2-4) implies d { F = (d t r)d r F + (diw)d u F + (^)^F , (^F) , (d^F) in the 
right-hand side of the above equations are given by d a F and d a dpFo (a, (3 = r,oj, fx). 
Therefore, our task is to find the concrete expressions of d a Fo and d a dpFo in terms 
of relevant geometric quantities. 

2.2. The CMB dipole formula 

First, we will obtain the CMB temperature dipole formula. For this purpose, 
we derive the expression for d a Fo. Our stating point is the Boltzmann equation for 
F (t,r,u,n), 

^F = d t F + fd r F + cod^Fo + (id^Fo = 0, (2-14) 

where r,uj,fi are defined for a given null geodesic curve 7. By differentiating this 
equation by a = r,u,fi, and using the formula (2-8), we obtain the first order 
differential equations for d a F, 



(2-15) 



The set of these equations can easily be integrated along the given photon trajectory 
7 to yield the solutions 

ud u F = M^o),, (2-16) 
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_ P -P(t,ti 



'HdrFo^ + iudMi [ dhe-r^Aih), (2-17) 
hi 

+{cjd u F )i £ dhe-^^ J B(h) - (j) J* dt 2 e~ p ^A(t 2 ) 1 , 



where 



and 



„ N N' S 



P(h,t 2 ) = T [ dt 



N' 



R S 



nY 

N f R' N'\ R S 



(2-18) 
(2-19) 

(2-20) 



« n f t2 , , N (R! N'\ R 
Wi'V^l dt \ ± ^{-R-N) + R-S 



i R 

ln iv 



1*2 



+ 2 



■Jti 



*2 



(2-21) 



and where the subscript 'i' denotes the value evaluated at the last scattering surface. 

Note that uid^Fo is constant for any null geodesies, as in (2-16). Furthermore, 
as shown in Appendix, by inspecting the regularity of F at the center, as well as 
the behavior of some relevant geometric quantities near the center, we can observe 
that d r Fo does not contribute to the leading behavior of diF in the limit r — > 0. 
Therefore, we only need to find the expression of d^F. 

Using the null geodesic equation for a radial geodesic, fi = ±1, dt/dr = ±S/N, 
we find that (2-20) and (2-21) become 



P(h,t 2 ) 

where 



In 



S_ 
N 



*2 



+ U(h,t 2 ), Q(h,t 2 ) = 2 



In 



R 

N 



*2 



+ 2U(h,t 2 ), (2-22) 



U(h,t 2 )= I dt 

It! 



N_S 
N~S 



(2-23) 



-2{7(t,*i)-{7(ti,ti) 



Substituting these into (2-18), we have 



(2-24) 
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Here we note that the second and third terms in the right-side of the above equation 
have the form 

where V(r,r\) is a function of r and r\ that is regular at < r < ri. For the radial 
geodesic that reaches r = at t = to, I(r) behaves as 

I(r) „ _I M (2 . 26) 

Substituting the above expression of d^Fo into 

diF~S ?-(f v ) = S ?-(^) , asr^O, (2-27) 



p p 
which is derived in Appendix (see (A- 16)), we have 



c2 i 

(OiF)o * Tirr 
N p 



\d r F )i + (ud u F )i { jT dre~ u ^A{t) - B 

(2-28) 

Thus, using this expression, we can write (2-12) as 



<5T\ (1) _ 5Ln-Q 



^e-uftoM ( JM) + P dre -u(t ,t) A{t) _ Bq 
Si \uduFoJi Jo 



(2-29) 

where 5Ln is the position vector of the observer. This is our dipole formula for the 
CMB temperature anisotropy in the most general spherically symmetric spacetime. 
The regularity of the metric, (2-1) at the symmetry center implies that A 2 ~ C\ + 
0(r 2 ), R/r ~ C 2 + 0(r 2 ), and S 2 ~ (R/r) 2 + 0(r 2 ) near the center, with d, C 2 
being some constants with respect to r. Using these estimations, we can check that 
the right-hand side of (2-29) is convergent, hence well-defined. 

2.3. The CMB quadrupole formula 

Next, we will derive the CMB quadrupole formula by inspecting the second- 
order derivatives, d a dpFo, (a, /3 = r,uj,/j,). By differentiating the first row of (2-15) 
with respect to lnw, r, and fx, we obtain 

±{(ojd^ 2 F } = 0, (2-30) 

^(wdudrFo) = T (J?) ood u d r F + A(ojdu) 2 Fo, (2-31) 

d , „ n „ ^ f N ( R' N'\ R S] ^ ^ 
«W = 2 { ± S ( R -n) + 5 -sj 

-juAAFo + B(wS„) 2 F„, (2-32) 
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for a radial geodesic for which /x = ±1. We can easily integrate this set of equations, 
and get the solutions 

(Lod w fF Q = {{udufFoji, (2-33) 

ududrFo = e- p ^{ud ul d r F ) t + {(ujd UJ ) 2 F } t [ dt ie ~ p(t ^ Aih), (2-34) 

hi 

ojd^Fo = e-^iud^Fo), - {ud w d r F } t J\t ie -^~ p ^^ (^j 

(2-35) 

Similarly, we can get the ordinary differential equations 
~(N\'-„ fN\" 



dt 



(d 2 r F ) = T 2 Hpj d 2 r F T (^) d r F + A'(ud u F )i + 2Aud u d r F , (2-36) 



|(<Wo) 



S \ R N J R S 
N\ 1 



drd^Fo - ^d 2 F 



+Bojd„d r F + AududpFo - (Jpj 8 r F + Cd^Fo + B'(ud u F )i, 

(2-37) 

1<^> = 4 { 4 (I " £) + 1 " 1} 8 2 F » " 4*Vb + ***** 

+D^F + 2 ^| - (ud u F )i, (2-38) 

for a radial geodesic, /x = ±1. We can also integrate this set of equations, and obtain 
the solutions 

d 2 F = e- 2P ^\d 2 F ) t 

+ jf <ttie- 2p ( t '* 1 ) |t 3 r F + A'(w0 w F o )i + 2Awd u d r F ^ , (2-39) 

ft ( AT 

+ J dt 1 e- p ^- Q ^ tl H--d 2 F + BududrFo + AudvdpFo 

- (f ) 5 r F + Cd M F + B'(wa WJ F )i | , (2-40) 

' ti 
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where 



+2Bud u d li F + DdpFo 



N (B! N' 



(2-41) 



ti 



^ = 2— — - — ± 6 - - - 



(2-42) 

Thus, we have the six solutions, (w<9 w ) 2 -Fo, ujd^drFo, etc. However, we note that 
(ujO^Fq is just constant. Furthermore, by inspecting the regularity of Fq at the 
center, as well as the behavior of some geometric quantities near the center, we can 
find that only 8 2 Fq becomes relevant to the evaluation of didjF in the limit r — > 0. 
In fact, as we show in Appendix (see (A-30)), 



didjF -> 2 U 



>i 3 



Si 



,p t p ] 



(f 2 )o + S 2 ^(d 2 r F ) 



+ 



p< 

W °ij + — T 

a± N J J p 2 



(udojFo)^, 



(2-43) 



where fi = d r iF (see (A-7) in Appendix), and a± = R/r, which corresponds to the 
'scale factor' perpendicular to the radial direction. Then, in terms of ji and 8 2 Fq, 
the CMB quadrupole formula (2-13) is written as 



5T\ (2) _ 1 Sx^xP 



2(5 lj - QiQj){h)o + ^^(9 r 2 F ) 



I "J. r , (nil II \ ^i^j 



S 2 



(ud u F )i 



+ 1 {(ST\ (1) \ {M-) 2i? o}, 



(wd u F )i ' 



(2-44) 



where J?j = SijX^ jr. So, the remaining task is to find the expressions of the leading 
behavior of f 2 and d 2 F$ at the center r — > 0. 

First, we note from (A-28) that (/2)o is given by 



(/2)o - 2 



1 /d r F 9„F f 



(2-45) 



In the limit r — > 0, from (2-18), the second term of this equation can be written as 



d„F v a ± 2 j\r 2 



iV 2 i? 2 



e -2U(to, ti )( diiFo) . T 



n N 2 e -2U(to,t) 

dr—r. n 



a ±0 N i ~2U(t Q ,tj) 

N 2 R 2 



(d^F )i ± 



a -LQ 

Jo r' a 

a ± 2 e- 2t/ (*°A) 
a± 2 



-d r F + -Bud^Fo 
S 



-8 r F + -Bujd u F 
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I i - 



-d r F + -Bcod^Fo 



± N 2 J ti r dt 



-2U(t ,t) 



-d r F + ^Bud u F 



From the second row of (2-15), we have 



d_ ( d r F 2U(t ,t) 
dt\ a 2 



S R N Na\ 



R N S a± 



A 



+ *-mto,t ){ujdu]Foh 



(2-46) 



N\'\ d r F 



2U(to,t) 



(2-47) 



Thus, (/2)o becomes 



a± 2 e -2l/(to,t0 



2r, : 



1 a , 2 f*o , AT2 



_a_L 
2 TV, 



o Jti 



dt- 



-8 r F + -BujO^Fq 



+ {- a - 2 [n-s)n b 



+ 



+a±- 



d f SB 



dt \Na 2 



(u)d u F )i 



~S 



d r F 



-2U(to,t) 



(2-48) 



Next, from (2-39), we obtain 
(d 2 F ) = e- 2P ^\d 2 F ) l 

+ jf ° dte~ 2P ^ S^AududrFo + A'(ud u F )i ± d r F j . (2-49) 

Thus, substituting (2-48) and (2-49) with (2-17) and (2-34) into (2-44), we finally 
obtain the quadrupole formula for the CMB temperature anisotropy. As in the dipole 
formula case, under the assumption that our metric (2-1) is regular at the symmetry 
center, we can check that the above quadrupole formula is well-defined. 



§3. The CMB temperature anisotropy in the LTB model 



In this section, the CMB temperature anisotropy formulae obtained in the pre- 
vious section will be given concrete expressions for the LTB cosmological models, 
and then will be applied in some specific LTB models considered in the numeri- 
cal analyses of Alnes and Amarzguioui, 53 ^ to check the consistency of the formulae. 
Further, the constraint on the location of off-center observers will be derived by our 
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analytic formulae with the latest WMAP data, thereby updating the previous results 
numerically obtained. But before doing so, we will first briefly recapitulate the LTB 
metric. 

3.1. The LTB spacetime 

A spherically symmetric spacetime with only non-relativistic matter, or dust, is 
described by the LTB metric, which may be given by setting in (2-1), 

N 2 = l, S= R 'f r ) , (3-1) 
1 — k{r)r z 

with k(r) begin an arbitrary function of r. Then, the Einstein equations reduce to 
f k\ 2 2GM(r) k(r)r 2 A . . M'ir) . . 

where Mir) is an arbitrary function of only r, and pit, r) is the energy density of the 
dust fluid. The solutions to (3-2) depend on the sign of k(r) and can be expressed 
in parametric form: For k(r) > 0, we have 

Mir) Mir) 
K\r)r {k(r)r 2 }2 

where t s (r) is an arbitrary function of only r. For k(r) = 0, we have 

R(t,r)= Q) 3 M§(r){t-* s (r)}i. (3-4) 

For k(r) < 0, we have 

Mir) Mir) 

R(t,r)= V / 2 (cosh7?-l), t-t s (r) = - J ^-'i S mh7 ] -r ] ). (3-5) 

— K\r)r |_/ c ( r ) r 2} 2 

The area radius R(t,r) vanishes at t = t s (r), so that t s (r) is called the big-bang 
time. The solutions admit three arbitrary functions k(r), Mir) and t s (r), but due 
to one degree of freedom in rescaling r, only two of them are independent. By 
appropriately choosing the profile of these two arbitrary functions, one can construct 
LTB cosmological models that can reproduce the observed SNIa distance-redshift 
relation. 

3.2. The analytic formula for the CMB dipole 

In the LTB cosmological model, in addition to a± = R/r defined previously, 
we also introduce the 'scale factor along the radial direction' by a//(t,r) = R'(t,r). 
Accordingly, we also define two Hubble expansion rates in the radial and azimuthal 
direction, respectively, by 



(3-6) 
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Prom (2-19) and (2-23), we obtain 

rt 2 



U(h,t 2 ) = - [ 2 dtH ih A = H' /n B = ±2(H /f - H ± ). 



(3-7) 



Then, the analytic formula for the CMB temperature anisotropy dipole (2-29) takes 
the form 



dT \m ( e -u(t , k ) 

— = TSLn-Ql 

T j | Si V w5 w Fo 



- ^ drH' u e- u ^ \ . (3-8) 



Now, using this formula, we derive some constraints concerning the position of 
off-center observers. In general, the CMB temperature anisotropy can be decom- 
posed in terms of the spherical harmonics Yi m by 5T/T = Yl<i m a i™Yim- We are 
interested in aio as the dipole moment. From (3-8), we obtain 



aio = T 



in 8L J e ~ U[t ° M) 
3" | Si 



-U(to,t) 



(3-9) 



Assuming that the universe is locally in thermal equilibrium at the last scattering 
surface, the first term in the bracket is of the order of {d r T/T)i because Fq is isotropic 
and depends only on w/Tj. In the models we consider in the present paper, this term 
can be neglected because the void size is sufficiently smaller than the horizon size 
and therefore, the observed region of the LTB universe is homogeneous with good 
accuracy on the last scattering surface. Furthermore, in order to estimate this term 
correctly, we have to specify the universe model before the last scattering. This is 
beyond the scope of this paper. Therefore, we have only estimated the contribution 
of the second term numerically. For the LTB model considered in 53), we have 
found that the induced aio is about 1.23 x 10~ 3 or less, according to WMAP data, 54 ) 
which implies that the distance from the observer to the center, 8L, has to be, 
5L < 16Mpc. This is consistent with the numerical result of 53). We have also 
applied this formula to various LTB models, and found, for example, 5L < 14Mpc 
in the Garfinkle model, 35 ) and 5L < 12Mpc in the GBH model. 39 ) 

3.3. The analytic formula for the CMB quadrupole 

As for the quadrupole moment in the LTB model, from (2-48) and (2-49), we 



obtain 



(/ 2 )o 



T iti e -mto,u) {dMi 



a 2 e -2U(t ,U) 

a±1 2ri 



-d r F + -jfBuduFo 



On 



to ^ 
dt- 



2( H// -H ±T ^) + (-t-)')<) r F {] 

a // a±J \a /n 



+ {-H'+2H // ^LB 



e 

d f SB 



dt \ Na 2 , 



{ujd w FQ)i 



-2U(t ,t) 



(3-10) 
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a 2 Fo)o = e^^(9 r 2 F ) J 

+ jf ° dte~ 2P ^ 1 2H' // ujd UJ d r F Q + H;j(ud u F )i + (J-^j d r^ 



(3-11) 



where ao = So = O// = a _i_0) an d £ = yl ~~ k(r)r 2 , just for notational simplicity, 
and where 



P(t 1 ,t 2 ) = T £ 2 dt(J-^J , (3-12) 



which is obtained from (2-20). Thus, we now have the analytic formula for the 
quadrupole moment of the CMB temperature anisotropy in the LTB model: (2-44) 
together with (3-10) and (3-11). 
From (2-44), we derive 



16vr {5Lf f 2(/ 2 ) , (d?F ) , (5" - <)„ 



, (aio) 2 (W) 2 fo), 

2^ (wa w F )i • 1 ; 

If the universe is locally in thermal equilibrium at the beginning, we can set {d a d^Fo)i 
to be zero. Further, for the same reason as we explained for the dipole formula, we 
neglect the term (d a d r Fo)i in the present paper. Under these assumption, we have 
estimated the quadrupole moment using this formula numerically for the model in 
53), and found that a 2 o — 8.61 x 1CT 7 . This is consistent with the numerical result 
of 53). For other models, for example, a 2 o — 5.51 x 10~ 6 in the Garfinkle model, 35 ) 
and a 20 ^ -9.27 x 1(T 7 in the GBH model. 39 ) 



§4. Summary and discussions 



In this paper, we have derived the analytic formulae for the dipole (2-29) and 
quadrupole (2-44) moments of the CMB temperature anisotropy in general spheri- 
cally symmetric spacetimes, including the LTB cosmological model as a special case. 
The formulae can be used to compare consequences of the LTB/local void models 
with observations of the CMB temperature anisotropy rigorously. The formulae also 
enable us to identify physical origins of the CMB temperature anisotropy in the LTB 
models. For example, in the CMB dipole formula (2-29), the first term comes from 
the initial (spherical) inhomogeneity at the last scattering surface, while the sec- 
ond term represents the integrated Sachs- Wolfe effect. Note that the first term also 
contains a contribution that reduces to the second-order ISW effect in the spatially 
homogeneous case. 

We have checked the consistency of our formulae for both dipole and quadrupole, 
with the widely-used recent numerical results for the special LTB model by H. Alnes 
and M. Amarzguioui. 53 ) Furthermore, we applied our formulae to other LTB models, 
such as those in 35), 39) and in particular, for the dipole moment, we found the 
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constraints on the distance between the void center and an off-center observer, by 
using the latest WMAP data. 

We can also utilize our analytic quadrupole formula to discuss the relevance of 
the LTB model to observed anomalies. For example, the observed magnitude of the 
quadrupole is known to be significantly lower than the ylCDM model predicts. This 
is usually understood as a cosmic variance, i.e., to be produced by a special feature 
of our Universe, one particular realization of the statistical ensemble. Because the 
local void model is one of such realization with a very low probability in the standard 
ylCDM model, it is tempting to see whether the quadrupole anomaly of the CMB 
anisotropy can be explained by a local void model. Unfortunately, however, the above 
analysis of the constraint on the observer offset by the dipole moment implies that the 
observed anomaly cannot be explained solely by the induced quadrupole moment in 
LTB models such as those in 35), 39), 53). Nevertheless, this result is not conclusive. 
For example, we have implicitly assumed that the off-center observer stays at a fixed 
comoving position. If the observer has a peculiar velocity pointed toward the center 
of the void, however, the value of 5L could be chosen to be much larger than the 
case with no peculiar velocity. If it is the case, then the observed anomaly of the 
quadrupole could be explained within the LTB models of 35), 39), 53). Therefore, it 
would also be worth attempting to develop other analytic formulae concerning CMB 
polarizations, lensing effects, etc. (cf. 55)) that can be used to distinguish the LTB 
and FLRW cosmologies. 
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Appendix A 

The regularity and derivatives of F$ near the center 

In this appendix, we discuss the regularity of the distribution function F at 
the symmetry center, and find which of d a F (resp. dadpF), a,f3 = r, become 
relevant in the leading behavior of the first- (resp. second-) order derivatives of F 
in the limit r — > 0. 

Mathematically, the distribution function can become singular at some radius 
including at the center. Furthermore, some authors concluded that a C 2_ -class 
singularity of the metric should be allowed at the center in order to construct a 
model with accelerated expansion exactly at the center. 36 ) However, such a model 
can be easily made smooth by appropriate smoothing and the original singular model 
can be recovered as a limit such that the smoothing length approaches zero. Hence, 
the final formulae for the dipole and quadrupole anisotropy of CMB given in the 
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present paper can be applied also to models with singularity in the metric or the 
distribution function if the results are finite. Hence, in this appendix, we assume 
that the metric and the distribution function are regular and smooth everywhere in 
Cartesian-type space coordinates on which the rotational symmetry group acts as 
on the standard Cartesian coordinates of the Euclidean space. 

In such a coordinate system x = (V), the spatial part of the metric (2-1), 
denoted here by g^, can be written as 

gij = S 2 Qi Qj + a\{5ij- fiiOj), ( A • 1 ) 

where Q l = x % jr and a± = R/r. For any smooth and rotationally invariant function 
h in this coordinate system, when it is Taylor expanded as 

h = ho + h i x i + h ij x i x j H , (A-2) 

each coefficient hij... must be a rotationally invariant constant tensor. Hence, these 
coefficient tensors vanish for odd ranks and can be written as the sum of products 
of the Kronecker delta Sij for even ranks. This implies that h is a smooth function 
of r 2 = 5ijX l x^ . 

Similarly, a smooth distribution function in this coordinate system can be Taylor 
expanded as 

F(t,x,p) = b (t,p) + b i (t,p)x i + b ij (t,p)x i x j + ■ ■ ■ , (A-3) 

where p = (p l ). When F is rotationally invariant, each term on the right-hand side 
is rotationally invariant separately. This implies that each b^...^ is an SO (3) tensor 
depending only on the non-trivial vector p l and therefore, can be written as the sum 
of products of the Kronecker delta 5^ and the vector p = (p l ). This implies that F 
can be written 

F(t,x,p) = /(t,x i x i ,|p|,x i p < ), (A-4) 
where XiX 1 = r 2 and \p\ = (5ijp l p^) 1 ^ 2 . Here, from 

P 2 = QijPY = Cixip 1 ) 2 + alSijpY, (A-5) 

where C is a smooth function defined by 

S" 2 - a 2 

C ^ r ) = —^^ ( A - 6 ) 

it follows that \p\ is a smooth function of t, oj = p, xip 1 and r 2 . Therefore, the 
regularity of F at the center is equivalent to the condition that the corresponding 
function Fo(t,r,u, p) can be written in terms of a smooth function / with the four 
arguments t,r 2 ,y = lnco,u = rp = Sxip 1 jp 

F (t,r,u,n) = f(t,r 2 ,y,u). (A-7) 



Then, 



8 t F = dif = {d l {r 2 )}f 2 + (d iy )f y + {d iV )f v , 
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where fa = d r 2f, f y = d y f, and f v = d u f . Here, because the spatial derivative of p 2 
can be written 



AY 



^-{rp r f + WLt-6 jk]f > p '- > ,•' + 2( '(>,/>' 



we have 



d t ( P 2 ) - 

d t (r 2 ) = 2x\ 

y riV ra± \2r ra± J \ b / 

dl u = S P - + rS' X - P --rSf d -^. 

2p a 



x l + C 



fir p 
S p 



p r p 

In particular, in the limit r — > 0, we see 

C -> {a ± (5" - a'i)} , 

^(p 2 ), ^(r 2 ), ^(ln W )^0(r), 



So P - 



P 



(A-9) 

(A-10) 
(A-ll) 

(A-12) 



(A- 13) 
(A- 14) 

(A-15) 



Therefore, from (A-8), we find that the first derivative of Fq behaves at the center 
as 



d lF ^ sAu)o = s p -( 9 » Fo 
p 



(A-16) 



P \ r Jo 

Next we study the second order derivatives of F with respect to x l , which are 
written as 

didjF = {^(r 2 )}/ 2 + {d,{r 2 )} [{d i (r 2 )}fa 2 + {d lV )fa y + (d^)f 2l/ ] 
+(didjy)f y + (djy) [{di(r 2 )}f y2 + (diy)f yy + ($!/)/„„] 
+{didjv)U + (^) [{^(r 2 )}/^ + (^y)/^ + ($!/)/„„] . (A-17) 

We find that 



did 3 {r 2 ) = 2<%, 



(A-18) 



y ) (V) 2 + 



a 



x l h 2—{rp r ){p l x : > + p>x l ) 



+ {-(r P n 2 + {a2±y 



didj(\noj) 



5 kl p*p l \5 t3 +2Cp l p>, 

r r 
2\ a (J2\a (Jl 



(A- 19) 



2p 2 



2(p 



,2\2 



N'Y x 
rN 



-x J 



rN 



Si 



(A-20) 



didjv 



r 



x l p 3 + 2x°p % rp 
p 2p 



+ x*d 3 (p 2 )} 



5 



+ J»W)} - S^r { d^ip 2 ) - -L- dl ( P 2 )d 3 (p 2 ) 



2p 



2p 



2p 2 



(A-21) 
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In particular, in the limit r — > 0, 

didjir 2 ) ^2<%, (A-22) 

''HV) "> 2 ( : ~- ) P% + 2C pV . (A-2:V) 

pV 7 





) P 2 5 







<_ 









(A-24) 



^^•(i/) -> 0(r). (A- 25) 

From these we find that in the limit r — > 0, 

a^-F -> 2%(/ 2 ) + { ( J - ^) <% + Co^} (f y )o + S 2 ^(U) . (A-26) 

Now, from (A-7), we find that 9 r Fo = 2r/ 2 + m/^j 9^Fq = rf u , ujO^Fq = f y , and 
8 2 Fq = 2/2 + 4r 2 /22 + 4r^/2^ + l~i 2 fw, and hence in the limit r — > 0, 

(A-27) 











1 /a r F 


2 V r 


fvv -> 













(A-28) 



(A-29) 



Thus, we finally obtain 



1 

2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 



didjF 2 r«5 ij - 5 2 ^) (/ 2 ) + S 2 ^(d 2 F ) 

+ { " ^r) ^ + Co ^} ( ^ Fo) - (A ' 30) 
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